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Abstract 

Composition fluctuations in disordered melts of symmetric diblock copolymers are studied by Monte 
Carlo simulation over a range of chain lengths and interaction strengths. Results are used to test three 
theories: (1) the random phase approximation (RPA), (2) the Fredrickson-Helfand (FH) theory, which was 
designed to describe large fluctuations near an order-disorder transition (ODT), and (3) a more recent renor- 
malized one-loop (ROL) theory, which reduces to FH theory near the ODT, but which is found to be accurate 
over a much wider range of parameters. 
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Diblock copolymers are linear polymers that contain two blocks of chemically distinct 
monomers. Diblock copolymer melts can exhibit both a disordered phase and a variety of peri- 
odic ordered phases. 1 Composition fluctuations in the disordered phase can be measured by small 
angle x-ray and neutron scattering experiments. The structure factor S(q) measured in such exper- 
iments is well approximated for very long polymers far from any transition by the random phase 
approximation 2 (RPA). The RPA for S(q) is based on the same physical approximations as those 
underlying the self-consistent field (SCF) theory that is often used to describe the ordered phases. 
Several closely related coarse-grained theories that attempt to improve upon the RPA 313 yield 
corrections that decrease with increasing chain length N, and thus reduce to the RPA in the limit 
iV — > oo. The simulations presented here were designed to test the quantitative accuracy of these 
theories. 

We consider the structure factor 5(g) in a dense liquid of symmetric AB diblock copolymers, 
each of N monomers. Let 



where 5ip(r) = ^ c a(i") — 5cb(t) and q = |q|. 5cj(r) is the deviation of the number concentration 
Cj(r) of i monomers from its spatial average. In a diblock copolymer melt, S(q) is maximized at a 
nonzero wavenumber q*. 

The RPA predicts 2 an inverse correlation function of the form 



where R g o = (N/Q) l l 2 b is the radius of gyration of a random walk polymer with statistical segment 
length b, and c is the total monomer concentration. The parameter Xe is an effective interaction 
parameter that is used in the RPA and SCF theories to characterize the degree of incompatibility. 
Here, F(x) is a known dimensionless function 2 that has a minimum at a value x* , yielding a 
maximum in S (q) at a corresponding wavenumber ql = x* / R g o. The subscript of '0' is used to 
denote RPA predictions for the quantities S(q), R g and q*. For symmetric diblock copolymers, 
F(x*) = 21.99 and x* = 1.95, so that S (q) diverges at a predicted spinodal value (xeN) s = 
10.495. 

A series of closely-related coarse-grained theories have attempted to improve upon the RPA by 
taking into account fluctuation effects that it neglects. 313 Here, we compare simulation results to 
the Fredrickson-Helfand 3 (FH) theory, which was the first such theory for diblock copolymers, and 
the renormalized one-loop (ROL) theory,— ~— which is the most recent. These theories all yield 




(1) 



cNS^(q) = F(qR g0 )-2x e N, 



(2) 



2 



predictions for S~ 1 (q) as a sum S^ 1 (q) = So 1 (q) + 5S~ 1 (q), in which the correction 5S~ 1 (q) 
is proportional to a small parameter iV~ 1/2 with N = Nb 6 c 2 . Physically, iV 1//2 = A^ 1 / 2 6 3 c is a 
measure of overlap: It is proportional to the number N 1 / 2 ~ R 3 /V of chains of excluded volume 
V = N/c that can pack in the volume R 3 pervaded by a chain of size R ~ y/Nb. 

The FH theory for S(q) was designed to describe the dominant effects of strong composition 
fluctuations very near the order-disorder transition (ODT) of a melt of long, symmetric diblock 
copolymers. It yields a simple nonlinear equation for the peak intensity S(q*). This can also be 
expressed in terms of an "apparent" interaction parameter x* a mat we define by fitting S(q*) to the 
RPA, by taking 

cNS- 1 (q*) = 2{(xeN) s - X :N}, (3) 
where (% e iV) s = 10.495. The FH theory predicts 

^ g ) = ^ _j_ g , (4) 

with B = 280. 

More recently, several groups have developed a family of renormalized one-loop (ROL) theo- 
ries of correlations in polymer liquids^, with a potentially wider range of validity. The character- 
istic feature of these theories is the development of methods to distinguish non-universal effects of 
monomer-scale correlations from the universal effects of long-wavelength correlations, and to ab- 
sorb the effects of short- wavelength correlations into renormalized values of the RPA parameters 
b and % e . Beckrich, Wittmer and colleagues in Strasbourg 7 8 have developed a theory for universal 
deviations from random-walk statistics in homopolymer melts, and verified their predictions by 
extensive simulations. Wang 6 gave predictions for S(q) near q = in homopolymer blends. Our 
group 10-14 has developed predictions for S(q) and single-chain correlations in both blends and di- 
block copolymer melts for arbitrary q, which reduce to the results of other authors in appropriate 
limits. The ROL theory for symmetric diblock copolymer s 1Qil2iB predicts a correction of the form 

cN5S~\q) = N- l l 2 H(qR gQl x f a N) , (5) 

in which the dimensionless function H is defined a sum of Fourier integrals that we evaluate 
numerically. This expression has been shown^ 2 *^ to approach the FH prediction of Eq. © very 
near the ODT. This ROL prediction also appears 1012 to be the first correction to the RPA within a 
systematic expansion of cN5S~ 1 (q) in powers of N^ 1 / 2 . Because the validity of this expansion is 
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not restricted to the vicinity of the ODT, we expect the ROL theory (but not FH theory) to remain 
accurate for N ^> 1 even far from the ODT. 

The limitations of the RPA have been well documented by previous simulations of diblock 
copolymer melts 1518 and experiments.— - — Both experiments and simulations show a decrease in 
q* with decreasing temperature T and a nonlinear dependence of the peak intensity S(q*) on 1 /T 
near the ODT that are not predicted by the RPA. Previous simulations have not, however, provided 
very precise tests of the absolute accuracy of the RPA or (particularly) of theories that predict 
corrections to the RPA. One reason for this has been the absence of a clear prescription for relating 
the interaction parameter Xe that is required as an input to these coarse-grained theories to the 
more microscopic parameters that are controlled in a simulation. The analysis presented here uses 
several methods to avoid or minimize this ambiguity. 

Our simulations use a potential energy similar to that of Grest and coworkers. 17 Non- 
bonded beads interact via a purely repulsive Lennard- Jones potential, with Uij(r) = 
Eij [4(r/er)~ 12 — 4(r/a)~ 6 + 1] for r < r c , with r c = 2 1/6 o\ Consecutive beads within each chain 
interact via a harmonic bond potential k(t — lo) 2 /2. All simulations discussed here use parameters 
£aa = £bb = kT, l = a, and k = 400/cTV -2 , with a total concentration c = 0.7cr~ 3 . The 
magnitude eab of the AB repulsion is controlled by a parameter a = Eab — £aa, which indirectly 
controls Xe- 

We have simulated melts of chains of N =16, 32, 64, and 128 beads in a periodic L x L x L 
cubic box. Two values of L were used for each chain length N, to monitor finite size effects. 
For each iV and L, we conducted parallel replica-exchange MC simulations^ 2 - of systems with 
different values of a at constant T. To sample configurations, we combined a hybrid MD/MC 
move, 23 in which short constant energy MD simulations are used to generate proposed MC moves, 
with reptation and double-rebridging 24 configuration bias moves. 

To test the RPA, FH, and ROL theories, one needs an unambiguous way to determine values for 
the parameters b and Xe that these theories require as inputs. Our comparisons with theory all use 
a value of b defined, as by the Strasbourg group, 8 as a limit b 2 = liniTv^oo 6R 2 (N) / TV '. This has 
been evaluated by numerically extrapolating results for Rg(N) for homopolymer melts (a = 0) of 
varying N, giving b = 1.41a. This definition of b is required for consistency with the ROL theory, 
which predicts that random- walk statistics 7 and the RPA for S(q) 10 with renormalized values of 
b and Xe are exact only in the limit iV — > oo. An analogous procedure for estimating for Xe{ a )> 
discussed below, is also used in some comparisons. 
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FIG. 1: q dependence of S(q) for N = 64, at a = 0, for L = 30 and L = 43. The RPA and ROL 
predictions are shown by dashed and solid lines, respectively. 

Limit Xe = 0: We first consider the special case of a system with a = 0, or 6ab = £aa, 
corresponding to \e = 0. This is an idealization of a neutron scattering experiment in which A and 
B blocks are labelled by differential deuteration, but are otherwise identical. In this limit, both the 
RPA and ROL predictions for S(q) depend only on intramolecular correlation functions for chains 
in a dense melt. The RPA, which assumes random walk statistics, predicts cNSo l (q) = F(qR g0 ). 
The ROL theory yields predictions that differ from the RPA only as a result of predicted 0(N~ l l 2 ) 
deviations from random walk statistics.— 

Fig. Q] shows simulation results and predictions for S(q) vs. qR g o in systems with = 64 
(N = 240) and a = 0. Here, R 2 g0 = Nb 2 /6, with b = 1.41er. Equivalent results are obtained 
for box sizes L = 30 and L = 43. The RPA (dashed line) slightly underestimates both the peak 
amplitude S(q*) and peak wavenumber q*. The ROL prediction (solid line), however, fits the data 
almost perfectly, with no adjustable parameters. We find a similar level of agreement at a = for 
other chain lengths. 

Peak wavenumber: We next discuss the evolution of peak wavenumber q* with changes in a, or 
X e (a). Fig. [2]shows our results for q* /q$, where q^ = 1.95/R g0 is the RPA prediction. The ratio 
q* /<7q is plotted vs. the quantity x^N that is defined in Eq. §3§ by fitting the observed peak intensity 
to the RPA. This plot thus shows peak wavenumber plotted vs. a measure of peak intensity. 

Both q* and S(q*) were determined by fitting values of S(q) at a discrete set of allowed 
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FIG. 2: Peak wavenumber q*, nomialized by the RPA prediction q$, vs x* a N for N =16, 32, 64, and 128. 
Data for two values of box size L are shown for each chain length N. ROL predictions are solid lines. 

wavevectors to a smooth function. For N =64 and 128 and large values of a, this procedure 
becomes unreliable because there are too few allowed values of q within the peak. In Figs. [2] and 
we thus show results for q* and S(q*) for each iV and L only over the range of values of a for 
which the peak remains broad enough to allow a reliable fit. Within these regions, results from 
different box sizes are consistent, and we see no evidence of a transition to an ordered phase. 

Simulation results for q* decrease monotonically with increasing x« or a, exhibiting a total 
drop of 15-20 % over the range shown, in agreement with previous results. 18 Note that g*/<?o > 1 
for small x* a N, consistent with the behavior shown in Fig. 1 for a = 0, and that q* /q£ < 1 only 
for x* a N ~ 5. The RPA predicts a constant q* / ql = 1 (dashed line). ROL predictions (solid lines) 
were obtained from parametric plots of predictions for q* and S(q*) over a range of values of XeN. 
The ROL theory captures all of the observed qualitative features, and is strikingly accurate for the 
longest chains (N = 64 and 128). The predictions appear to become systematically more accurate 
with increasing N, consistent with the claim 1011 that the ROL theory is the first correction to the 
RPA within an expansion in powers of N^ 1 ^ 2 . 
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Estimating Xe(«)-' In all of the theories considered here, the degree of incompatibility between 
A and B monomers is characterized by an effective interaction parameter Xe- To test the predic- 
tions of these theories for how the peak intensity S(q*) evolves with increasing % e , we need an 
independent estimate of how Xe{ a ) depends on a. (The comparisons shown in Figs. 1 and 2 did 
not require this.) 

We estimate Xe(ct) using a method based on thermodynamic perturbation theory. Ref. 14 dis- 
cusses a perturbation theory for a structurally symmetric blend of A and B homopolymers, both 
of length N, with a pair potential Uij(r) = tiju(r), with tAA = ^bb and a = cab — ^bb, in which 
the excess free energy per monomer f ex (a,N) is expanded in powers of a. In such a model, 
fex{& — 0, N) = 0. It was shown 14 that f ex is given to 0(a) by f ex (a, N) ~ az(N)(j) A (j) B , where 
(pi is the fraction of i monomers, and where 

z(N) = J dr g(r;N)u(r). (6) 

Here, g(v; N) is the inter-molecular radial distribution function in the reference state with a = 0. 
By comparing this perturbation theory to ROL predictions for S(q) at q = 0, which reduce to the 
RPA in the limit iV — > oo, while also expanding Xe(a) in powers of a, it was found that the RPA 
parameter x e is given to O(a) by 

Xe - z(oo)a/kT, (7) 

where z(oo) = limA^oo z(N). For the model considered here, 14 we find 2(00) = 0.2965. Pertur- 
bation theory provides a useful description of the homogeneous state of both blends and diblock 
copolymer melts because the critical values of a and Xe are proportional to 1/N, implying that the 
accuracy of perturbation theory must improve with increasing N. 

Peak intensity: Fig. |3] shows how the inverse peak intensity cNS^ 1 ^*) varies with Xe-^> 
plotted using Eq. © for Xe(cn). Fig. |3(a)| compares results for chains with iV = 64 to the RPA 



(straight dashed line), FH (dot-dashed), and ROL (solid) predictions. Fig. |3(b)| compares results 
for chains of length iV = 32, 64, and 128 to RPA and ROL predictions. The comparisons to theory 
involve no adjustable parameters. 

Several qualitative features are immediately apparent: The quantity cNS -1 ^*) is greater than 
the RPA prediction (suppressed fluctuations) for x e iV > 6, as predicted by the FH theory, but 
is less than the RPA prediction (enhanced fluctuations) for XeN ^ 6. For XeN ^ 10, results 
for different converge towards the RPA prediction with increasing N, as predicted by both FH 
and ROL theories. In Fig. |3(a)[ the FH and ROL predictions for N = 64 are very similar to 
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FIG. 3: Inverse peak intensity vj. XeN for = 64 (panel a, upper) and for iV = 32, 64, and 128 (panel 
b, lower). For each N, results are show for two values of the box size L. Xe is estimated using Eq. (0, 
using z c = 0.2965. Straight dashed lines show the RPA prediction. Solid lines are ROL predictions. The 
dot-dashed curve in panel a is the FH prediction. 



each other and to the simulation results in the strong fluctuation regime, where XeN ^ 15 and 
ciVS'~ 1 (g*) ^ 1. This is the regime that the FH theory was designed to describe. For XeN ^ 15, 
however, the FH theory fails, while the ROL theory remains remarkably accurate, for this chain 
length, down to XeN = 0. The ROL theory correctly predicts a change in the sign of the deviation 
of cNS~ 1 (q*) from the RPA prediction at a value of XeN — 6, which is shown in Fig. |3(b)| to be 
almost independent of N. The FH theory incorrectly predicts a rather large positive deviation for 
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all XeN > 0. In |3(b)[ the ROL theory is seen to be quite accurate for all three values of N for 
X* a N ;$ 7, and to remain accurate for larger values XeN for N =64 and 128. ROL predictions 
differ significantly from the results only for the shortest chains shown, with iV =32, at large values 
of XeN. Results for N = 16 are not shown, but exhibit the same trends, and differ even more from 
the ROL prediction. The apparent tendency of the ROL theory to become more accurate with 
increasing iV is consistent with the claim that it is part of a systematic expansion in powers of 

JV-V2. 

The comparison shown in Fig. 3 is subject to at least two types of error: (i) Errors in the ROL 
theory itself, due to truncation of the loop expansion at first order in iV -1 / 2 , and (ii) Errors in our 
estimate of % e , due to the truncation of the Taylor expansion of Xe(«) at first order in a. The 
comparisons shown in Figs. 1 and 2 are subject only to type (i). The error in Xe{ a ) ls expected 
to depend only on a, and to increase with a. It should be most serious for large values of XeN 
and small N, because the value of a required to obtain a specified value of XeN increases with 
decreasing N. The fact that the largest discrepancies between the ROL theory and the simulation 
results occurs in Fig. 3 for the shortest chains at large values of XeN suggests that our use of a 
simple linear approximation for Xe( a ) mav be a significant source of error for the shortest chains. 

In summary, we have presented an unusually systematic simulation study of how S(q) depends 
upon both chain length and degree of AB repulsion in AB diblock copolymer melts, and com- 
pared our results to several theories. The use of theoretically motivated procedures for obtaining 
independent estimates of the RPA parameters b and Xe, defined by extrapolation to N = oo, al- 
lowed much more precise comparisons to theory than has previously been possible. The range of 
values of N studied here, N < 480, overlaps the lower end of the range studied in experiments 
on polymers, for which N is more often 10 3 — 10 4 . A recent theory, the ROL theory, is shown 
to be a substantial improvement over the FH theory, with a much wider range of validity. The 
ROL theory is strikingly accurate for the longest chain lengths studied, and for sufficiently small 
values of XeN. All our results appear to be consistent with the claim that this theory is the first 
correction to the RPA within a systematic expansion, and should thus become more accurate with 
increasing N. The results indicate that theories developed to describe subtle deviations from the 
random-walk and RPA theories in liquids of very long polymers, based on an expansion in powers 
of iV -1 / 2 , can accurately describe liquids of surprisingly short chains. 

This work was supported by NSF grant DMR-097338. 
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